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Simulations de l'écoulement turbulent marin avec un 
modèle de deconvolution 

Anne-Claire Bennis[ï], Roger Lewandowski H et Edriss S. TitiH 

. Résumé 

O 

On introduit une équation de déconvolution qui généralise l'algorithme de Van-Cittert 
£Li pour des conditions aux limites de type océananique avec vent fixé. On en déduit un modèle 

de SGE pour lequel on a existence et unicité d'une solution régulière. Nous détaillons un 
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ensemble de simulations numériques qui montrent l'intérêt pratique du modèle. 

Abstract 

We display a continous équation for the déconvolution process that generalizes the 
Van Cittert algorithm in the case of oceanic boundary conditions for a given fixed wind. 
We deduce a LES model for which we have existence and uniqueness of a strong solution. 
Finally, we display several numerical simulations showing the practical interest of the model. 



Abridged English version 

A LES model based for the Navier-Stokes équations (NSE), where the convection term is 
l/"} ' replaced by (H^(u) ■ V)u, was introduced recently in [5]. In this model, H^(u) dépends 

on u in the following way. For a given u, with V • u = 0, we filter u by solving the Stokes 
(Helmholtz) like elliptic problem — a 2 Au+u+V7r = u, Vu = subject to the appropriate 
boundary conditions, and we then set Au = u. We then follow the Van Cittert 's algorithm 
Q ■ (see [5]), uo = û", u n+ i = u n + {û — A _1 u n } and we set u^r = i?jv(u) for some N. This 

model makes sensé in the case of periodic boundary conditions (see [S]). Nevertheless, 
we cannot use it in the case of realistic nonhomogeneous boundary conditions, such as 
those satisfied by oceanic flows with a fixed wind stress at the surface under the rigid lid 
^ . hypothesis. For this reason a substantial modification is necessary. 

Our contribution starts from observing that the Van Cittert algorithm is in fact nothing 
but a finite différence scheme of a certain évolution équation. Taking r as the phony time 

évolution parameter, this équation for scalar quantifies takes the form —a 2 A— — h v = u, 

_ OT 

v\ t= q = u. In this work, we make use of this remark to study the case of an oceanic flow 
satisfying the boundary conditions (|2.ip . below. Here x = (x,y,z) G IR 3 , O C H 2 is a 
bounded smooth domain, Q = Ox[-h,0], T s = Ox{0}, T t = dOx[-h,0], T b = Ox{-h}; 
r = r s ur;Urb. Finally t € [0, T] is the real time parameter. We first introduce the System 
(I2.4p - (l2.5j) . This yields the déconvolution model (12.7p . In this model, we do not make the 
hydrostatic balance assumption in order to capture and model the vertical convection. 
For the sake of simplicity, we do not consider the Coriolis force, which can be added to 
(12.7p . without causing any difficulties in our analysis. We now give a short description 
of our results that will be reported in détail in pp. The space functions are defined in 
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(|2.2p below. For a given wind stress V, -0 = tfjÇV) is a divergence free vector field which 
satisfies the boundary conditions (|2.1|) and is defined by (|2.3p . Let r > be given, and 
u = u(i, x) with V • u = 0. We set H T (u)(t, x) = v(r, t,x) +i/;(V)(x) where v(r, t,x) is 
solution of System (|2.4|) - (|2.5p for a fîxed time t (see (j2.6|) below). Notice that when r = 0, 
-ffo( u ) = û and therefore our model is precisely the Leray-a model, as considered in [2] 
and [3]. When r = iV is an integer, this model remains the same as the one introduced 
in [5] for the periodic boundary conditions case, and here it is adjusted to satisfy physical 
boundary conditions. Our main results are 

Theorem 0.1 Assume V G H 3 (0) 2 with a compact support, uo G IHi (see définition 
(\2.2\ ) and f G C°([0, T], IHi). Let t > be given, Problem (2.7) fias a unique strong 



solution (u T ,p T ) G L 2 ([0, T], H 2 {Vif) x L°°([0, T], fl" 1 ^)) with u T G C([0, T], L 2 (nf), 
d t u T G L 2 ([0,T] x fi) 3 . 

Theorem 0.2 There exists (r n ) ne 7V; which goes to infinity, when n goes to infinity, 
and such that the séquence (u Tn ) ng N converges weakly in L 2 ([0, T], IHi) and strongly in 
L 2 ([0,T],m ), to a field u = (u h ,w) G L 2 ([0, T], Mi) n L°°([0, T], L 2 {Qf), a Leray-Hopf 
weak solution to the NSE, with dtu G Ls ([0, T], IH^) and suc/i i/iai for ail v = (vh,6) G 
L 4 ([0,T],IHi) ([Oj) is satisfied. Moreover, u(i, •) converges weakly to uo in L 2 (f2) 3 as t 
goes 0. Moreover, u satisfies the energy inequality. 

Finally, we display numerical tests by using FreeFem++ [3] in the 2D case. Our code is 
semi-implicit and makes use of a method of the caracteristics adjusted to the deconvolution. 
Two cases have been considered. For each case, the wind stress is equal to V(x) = 
a(x) sin(-7rx), x = (x,y), x^ = x G T s ?» [0, 1], where a(x) has a compact support and it 
localizes the stress. In the first case, is a square, At = 0.2 s, a = 0.1m, T = 90 s. The 
DNS exhibits a good behavior. We compare the cases where the order of deconvolution r 
takes for values r = (Leray-a model), r = 5 and r = 20, At = 1. Vertical profils are 
displayed in Figure 13.11 They confirm the convergence Theorem and the validity of the 
code. The second case is the one where il has a topography (see Figure [372]) . At = 0.02 s, 
q = 0.1m, T = 10 s, r = 5, At = 1. Poiseuille Flows have been imposed on the latéral 
boundaries to enhance turbulence and the DNS converges very slowly. Figure 3.3 displays 
one profil and underlines the improvement of the computation using the model compared 
with the Leray-a model. Therefore, our model appears to produce good results and hence 
could be suggested as a good alternative for simulations in real situations. 



I Introduction 

On a introduit dans [5] un modèle de SGE basé sur les équations de Navier-Stokes, dans 
lequel le terme de convection vaut (Hn(u) ■ V)u, Hn(u) étant déduit de la vitesse u 
(avec V • u = 0) en partant de la filtration —a 2 Au + u + Vit = u, V • u = 0. On pose 
Au = u, que l'on déconvole à l'aide de l'algorithme de Van Cittert, uo = u, u n+ i = 
u„ + {u — A _1 u„},Uiv = Hn{u). Ceci a un sens dans le cas de conditions aux limites 
périodiques et permet d'approcher les équations de Navier-Stokes par des problèmes bien 
posés avec un bon contrôle de l'erreur quand les lois de Kolmogorov sont satisfaites. 

II semble difficile d'appliquer cette procédure telle quelle dans le cas de conditions aux 

limites réalistes, comme par exemple pour l'Océan avec une tension du vent en surface 

fixée, sous l'hypothèse du toît rigide. La remarque principale dans ce travail est que 

l'algorithme de Van Cittert peut être vu comme l'équation aux différences finies d'une 

du 

équation d'évolution qui s'écrit pour des scalaires sous la forme —a 2 A— — \-v = u, u\ T= o = 



n, où r est le paramètre de déconvolution. Nous exploitons cette remarque de manière 
rigoureuse pour les conditions aux limites (|2.ip satisfaites par la vitesse de l'écoulement 
marin, en introduisant le système (|2.4p - (|2.5p . Cela nous permet d'introduire le modèle de 
SGE (|2.7p . Dans ce modèle, nous ne faisons pas l'hypothèse hydrostatique pour mieux 
décrire les effets convectifs et capter la couche de mélange. Pour simplifier la présentation, 
nous n'avons pas pris en compte la force de Coriolis, qui peut être rajoutée dans (j2.7|) sans 
changer la nature de notre analyse. 

Nous montrons que ce modèle conduit à un problème bien posé et nous avons un résultat 
de convergence vers les équations de Navier-Stokes quand le paramètre de déconvolution 
tends vers l'infini. Dans une série de tests numériques, nous validons le modèle dans le cas 
2D en comparant les résultats obtenus avec ceux d'une DNS convergente. Nous produisons 
des calculs numériques dans un cas avec topographie, cas pour lequel la DNS a un mauvais 
comportement. Le modèle donne des bons résultats numériques pour r = 5 qui affinent 
ceux donnés par les Leray-a modèles habituels adaptés à nos conditions aux limites avec 
notre méthode, puisqu'ils correspondent au cas r = (voir les modèles d'origine dans [2], 
[3]). Cette note annonce des résultats développés dans l'article en cours de finalisation pQ. 



2 Résultats théoriques 

2.1 Géométrie et espaces fonctionnels 

On note x = (x,y,z) un point de M 3 , O C IR 2 un domaine borné régulier de M 2 , $7 = 
O x [—h, 0], où h > 0. Soient T s la surface, r s = 0x {0}, et T; et Tj, les frontières latérales 
et le fond, Ti = dO x [-h, 0], T b = O x {-h}. On note enfin r = T s U T t U r 6 . 
Soit u = (u/j, w) la vitesse de l'eau de mer, = (u, v) sa partie horizontale. Les conditions 
aux limites classiques sont (sous l'hypothèse du toît rigide) 

(2.1) w\ r = 0, u/j| ri ur 6 =0, Vu h ■ n\ Fs = V, 

où n est la normale sortante, V = V(x, y) est le vent (ici fixé) défini sur O. On identifie 

et T s . Les espaces fonctionnels sont les suivants : 

(2 2) Ho = {u = (u h , w) G L W; u • n| r = 0, V • u = 0}, 

1 • > IHx = {u = (u h , w) G W, Mv = 0, u h |r iUr() = 0, V • u = 0}, 

Nous montrons dans pQ l'existence d'une base spéciale adaptée, résumé dans ce qui suit. 

Lemme 2.1 77 existe une suite de réels strictement positifs (A n ) ng ]N des suites (e n )n£ti, 
(7r n ) ng N où e n S IHini/ 2 (r2) 3 , 7r n S H 1 ^) pour chaque n et tels que — Ae n +V7r n = À n e n . 

La démonstration utilise un argument de symétrie pour se ramener un problème de Stokes 
sur un cylindre avec des conditions aux limites de Dirichlet homogènes. Un résultat ana- 
logue avec des équations primitives a déjà été prouvé dans [7] avec la même technique. 
Nous montrons dans [T] que l'injection IHi C IHo est dense. Nous terminons cette section 
par le relèvement de la condition aux limites. On suppose que V = V(x, y) est à support 
compact dans T s . Soit rj) = ij)(x, y, z) le champ de vecteurs défini dans f2 par 

(2.3) %l>{x,y,z)=\ = =^(V), 

\tpv J \ -K{z){V h -V){x,y) J 

où k(z) = J\ p{z)dz et p{z) = j^z 2 + z + ^. Il est facile de vérifier que ip = tp(V) est à 
divergence nulle dans Q et satisfait les conditions aux limites (|2.ip . 



2.2 La Filtration 

On suppose que V G H 3 (0) 2 et est à support compact, a > est fixé, u G IH4. Soit 
J -a 2 Aïï + ïï + Vr = u - -0(V), V • H = 0, 
I ûJ|r = 0, û/i|r,ur b = 0, VÛ^ • n|r s = 0. 



(2.4) 



En raisonnant par symétrisation on montre dans [T] que (|2.4p admet une unique solution 
(ïî, r) avec u G i/ 3_e (f2) 3 pour tout e > 0. Soit û" défini par u = u + ifiÇV). On note 
que u a une trace sur T s ainsi que Vû" et donc u satisfait (|2.ip . Soit u = u(t, x) à 
divergence nulle. L'équation de déconvolution est, où r est le paramètre de déconvolution 
et v = v(t, t, x) = (v/i, 0), t G [0,T] est fixé, 

-a 2 A^+v + Vvr = u-^(V) = û, V • v = 0, 

vh|r ; ur b = 0, Vv h -n|r s =0, 9\ dn = 0, v(0, x) = û + <0(V) = n. 

Soit r > 0, H T (u)(t, x) = v(r, t, x) + •0(V)(x). Avec la base du Lemme [27Tt on écrit pour 
uGL 2 ([0,T],IHi)nL 00 ([0,T] ) IHo) ) u = E^i^(*)ejfe et </>(V) = £~ x ^e fc . On a 

(2.6) fl T (u)(t,x) = u(t,x)+y;(û fc (<)-u fc (t))ex P (~^-)e fc (x), n fc (t) = , 

fc=i 

2.3 Le modèle de déconvolution 

Le modèle de SGE que nous considérons est le suivant, avec u la vitesse et p la pression : 

| ftu + (£T T (u) • V) u - !/Au + Vp = f, V • u = 0, 
|w| r = 0, u h | riUr6 =0, Vu h -n| rs =V, u t=0 = u , 

Notons que quand r = 0, il s'agit du Leray-Q modèle ([2]). Nous prouvons dans [T] : 

Théorème 2.1 Soit V G iî 3 (C) 2 à support compact, u G Hi et f G C°([0, T], IHi). 
Powr r > /îxe, /e système l \2. 7| ) admet une unique solution (u T ,p T ) G L 2 ([0, T], H 2 (fl) 3 ) 
x L^QO^iï 1 ^)) avec u T = u r (t,x) G C([0, T], L 2 (fi) 3 ), d t u T G L 2 {[0,T] x fi) 3 . 

Théorème 2.2 II existe une suite (r n ) n£ Ar qui tend vers l'infini quand n tend vers l'infini, 
telle que (u Tn ) n€ ]N converge faiblement dans L 2 ([0, T], IH4), fortement dans L 2 ([Q, T], IH ), 
vers un champ u = (u^io) G L 2 ([0, T], Mi)nL°°([0, T], L 2 (fi) 3 ) avec d t u G Li([0, T], IHi) 
et tel que pour chaque v = (v h ,0) G L 4 ([0, T], IHi) on a 

(2.8) <d t u,v>- f [ u(g)u:Vv + i/ / f Vu:Vv- f [ V-v h = [ [ fv. 

Jo Jn Jo Jn Jo Jr s Jo Jn 

De plus, u(i, •) tend faiblement vers vlq dans L 2 (Q) 3 quand t tend vers 0, et vérifie 
l'inégalité d'énergie. 



3 Résultats numériques 

On résoud (12. 7p par la méthode des éléments finis en utilisant FreeFem++[l] dans le cas 2D, 
pour r fixé. On détermine la pression et la vitesse par un schéma semi-implicite initialisé 
avec des solutions d'un problème de Stokes. On a écrit un algorithme de résolution itératif 
qui adapte la méthode des caractéristiques au cas du transport déconvolé. On arrête la 
simulation lorsque les algorithmes ont convergé. On présente un premier cas servant à 
valider notre modèle de déconvolution et un deuxième cas qui montre son intérêt pratique. 



3.1 Turbulence induite par le vent de surface 



On considère $7 = [0, 1] x [—0.5,0], triangulé avec h max = 0.014 m et h m i n = 0.006 m. On 
choisit V(x) = a(x) sin(7rx) où a(x) est une fonction de localisation à support compact. On 
compare les résultats des simulations DNS et LES pour des paramètres de déconvoution 
successivement égaux à r = 0, r = 5 ou r = 20 grâce à des profils verticaux ; ici At = 1, 
a = 0.1 m, Ai = 0.2 s et on a fait 450 itérations en temps, ce qui nous donne le résultat pour 
T = 90 s. On trouve une erreur L2 moyenne de 0.5 % entre les résultats DNS et les résultats 
du modèle de déconvolution pour r = 20 et de 1% pour r = 5, ce qui est satisfaisant et 
ce qui valide notre modélisation. De plus, on a divisé l'erreur L2 moyenne d'un facteur 4 
entre les résultats du Leray-a modèle (r = 0) et ceux du modèle de déconvolution pour 
r = 20, ce qui montre l'amélioration apportée. 



Profils x=0.25m Profils x=0.75m Profils x=0.875m 
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Figure 3.1: Profils verticaux de la vitesse horizontale pour la DNS (DNS), pour le modèle de 
déconvolution avec r = 5 (LES-1) et avec r = 20 (LES-2) et pour le modèle Leray-a (Leray-a) 



3.2 Cas avec bathymétrie 

On considère le domaine O décrit ci-dessous avec une bosse, triangulé avec h max = 0.010 m 
et hmin = 0.003 m. 
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Figure 3.2: Géométrie du domaine à gauche et Résidus des équations de Navier-Stokes à droite 
(DNS en rouge, modèle de déconvolution en bleu et modèle Leray-a en noir) 



On prend le même vent à la surface que précédemment et on impose un flot de poiseuille 
sur les bords d'entrée et de sortie I\ et T; 2 afin d'augmenter le niveau de turbulence. On 
conserve les paramètres a = 0.1m, At = 0.02 s, r = 5 et on fait 500 itérations en temps, 
soit T = 10 s. On observe que la DNS converge difficilement avec des résidus élevés, ce 
qui n'est pas le cas des modèles LES (voir Figure 3.2). 

On observe que notre modèle améliore 
la solution donnée par le modèle Leray- 
a, ce qui est en accord avec les con- 
clusions précédentes. Ce profil mon- 
tre qu'il faut augmenter le nombre de 
déconvolution qui est égal à 5 comme 
préconisé par Adams et al [6]. Il est 
aussi nécessaire d'abaisser la valeur du 
paramètre de filtration a, afin d'avoir 
un temps de calcul raisonnable par un 
abaissement du taux de déconvolution. 
Ce cas montre que notre modèle peut 
corriger les effets de la filtration et mon- 
tre aussi que le calage des paramètres 
de déconvolution et de filtration n'est 
pas le même selon les cas étudiés. 
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Figure 3.3: Profil vertical de la vitesse horizontale. 
Modèle de déconvolution en bleu (LES-1), modèle 
Leray-a en noir (Leray-a) et DNS en rouge 
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